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ABSTRACT 

We present near-infrared H and i^-band spectra of the z = 6.41 quasar SDSS Jl 14816.64+525150.3. 
The spectrum reveals a broad Mgn A2799 emission line with a full-width half-maximium of 6000 kms -1 . 
From the peak wavelength of this emission line we obtain a more accurate redshift than is possible from 
the published optical spectrum and determine a redshift of z = 6.41 ± 0.01. If the true peak of the Ly a 
emission is at the same redshift, then a large fraction of the flux blueward of the peak is absorbed. The 
equivalent width of the Mgn emission line is similar to that of lower redshift quasars, suggesting that 
the UV continuum is not dominated by a beamed component. Making basic assumptions about the 
line-emitting gas we derive an estimate for the central black hole in this quasar of 3 x 10 9 solar masses. 
The very high luminosity of the quasar shows that it is accreting at the maximal allowable rate for a 
black hole of this mass adopting the Eddington limit criterion. 

Subject headings: galaxies: formation - quasars: emission lines - quasars: individual (SDSS 
J114816.64+525150.3) 



1. INTRODUCTION 

The vast energy requirements of the most luminous 
quasars can be met by a model invoking extraction of grav- 
itational potential energy from matter falling toward a su- 
permassive (~ a billion solar mass) black hole. The corre- 
lation between black-hole mass and galaxy mass observed 
in the local Universe strongly suggests that the most lu- 
minous quasars will reside in the most massive galaxies in 
the most massive dark matter halos. The luminous quasars 
at redshifts z > 6 found in the Sloan Digital Sky Survey 
(SDSS) therefore pinpoint the earliest massive objects to 
form (Fan et al. 2001, 2003). 

The recent discovery of SDSS J114816. 64+525150.3 
(hereafter SDSS J1148+5251) at a redshift of z = 6.43 
by Fan et al. (2003) makes it the most distant known 
quasar, observed only 840 million years after the begin- 
ning of the universe (cosmological parameters of Hq = 
70 km s™ 1 Mpc^ 1 , Qm = 0.3 and Q\ = 0.7 are assumed 
throughout). The very high luminosity of this quasar 
(M1450 = —27.8) suggests that it contains a black-hole 
mass of at least several billion solar masses. 

We have obtained near-infrared spectroscopy of SDSS 
Jl 148+5251 to search for the Mgn emission line (rest 
wavelength 2799 A) redshifted into the if-band. This line 
is particularly useful since it is usually found close to the 
systemic redshift, unlike higher ionization lines such as 
Civ which typically show large blueshifts ~ 1000 kms^ 1 
(Richards et al. 2002). The redshift for this quasar given 



by Fan et al. was based on the position of a heavily ab- 
sorbed Ly a line and the onset of the Lyman break and is 
therefore quite uncertain. An accurate redshift is impor- 
tant for determining the strength of the proximity effect 
and for molecular line searches. Additionally, it has re- 
cently been shown that the width of the Mgn emission 
line can be used to derive the mass of the quasars central 
black hole in the same way that the H (3 line has tradi- 
tionally been used (McLure & Jarvis 2002). 

We describe the spectroscopic observations and mea- 
surements made from them in §2. The importance of the 
new redshift determination is discussed in §3. We deter- 
mine an estimate for the mass of the black hole in §4. Our 
conclusions are presented in §5. 

2. NEAR-INFRARED SPECTROSCOPY 

2.1. Observations 

We observed SDSS Jl 148+5251 with the United King- 
dom Infra-Red Telescope (UKIRT) equipped with the 
UKIRT Imager Spectrometer (UIST) on 2003 January 10. 
UIST is a 1-5 /im imager-spectrometer with a 1024x1024 
InSb array. Conditions were photometric with seeing at 
-fCband approximately 0.6 arcsec. The HK grism was 
employed with a 0.48 arcsec wide slit giving continuous 
wavelength coverage from 1.4 to 2.5 microns with a resolv- 
ing power of ~ 500. The pixels have size 10.9 A spectrally 
and 0.12 arcsec spatially. The target was nodded along 
the slit a distance of 12 arcsec in between exposures to en- 
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Fig. 1. — Upper: UIST H and .K"-band spectrum of SDSS J1148+5251 (thick line; smoothed with a 7 pixel boxcar). The location of the 
broad Mgll emission line is marked. The smooth curve shows the best fit comprising a power law continuum (a = 0.2), a broadened Fen 
template and a broad Mgll line. Lower: The residual spectrum after subtraction of the power law continuum and Fell template (thick line). 
The best fit to the broad Mgll doublet (FWHM=6000 kms -1 ) is also shown here (thin line). The shaded region shows the ±1<t noise level 
of the observation. 



able accurate sky subtraction. The total integration time 
of the observation was 6480 seconds. 

The data frames were flat-fielded and then adjacent 
pairs were subtracted from each other to subtract the sky 
emission and produce positive and negative spectra of the 
target. These data were combined and any pixels devi- 
ating from the mean at the 3a level were excluded since 
this indicates contamination by cosmic rays. The resid- 
ual background level was subtracted from the combined 
image. The positive and negative quasar spectra were ex- 
tracted from this image using apertures of 1.0 arcsec and 
combined. 

Flux calibration and atmospheric extinction corrections 
were made by dividing the spectrum by that of a star 
of spectral type F6 which was observed immediately pre- 
ceding the quasar. An accurate wavelength calibration 
was made by comparison to an Argon arc lamp. This 
solution was then checked against the positions of about 
20 sky emission lines giving an rms uncertainty in the 
wavelength calibration of 1A. The absolute flux scale of 
the spectrum was set by comparison with the magnitude 
measured in a snapshot if-band image. The image gave 
K = 16.86 ± 0.05 consistent with the K' = 16.91 quoted 
by Fan et al. (2003). 

2.2. Emission line fitting 

The spectrum of SDSS J1148+5251 is shown in Fig. 
1. It shows a broad emission line at approximately the 
position expected for redshifted Mgll. To determine the 
wavelength and width of the Mgll emission line, a line- 
fitting process was applied which involves the determina- 
tion of the underlying continuum and Fe II emission sur- 
rounding the Mgll line. Due to the blending of the Mgll 
line emission and the Fe II features in the 2700 A - 2900 A 
region, the fitting of the continuum and Fen emission is 
confined to two bands either side of the Mgll line (rest- 
frame 2250 A - 2700 A and 2900 A - 3100 A). 

The Fe II emission in this spectral region is modeled us- 



ing an Fen template based on an archival HST FOS spec- 
trum of the narrow line Seyfert galaxy IZwl, which is 
notable for its extremely strong Fe n emission. Due to the 
fact that the Fe n emission in I Zw 1 is relatively narrow 
(FWHM ~ 900 kms -1 ), this template can be smoothed in 
order to match the much broader Fen features typical of 
powerful quasars. The fitting of the Fen emission tem- 
plate and AGN continuum are performed simultaneously, 
with the amplitude and FWHM of the Fe n iron template 
left as free parameters. The AGN continuum is modeled 
as a power law which is normalized with the continuum 
magnitude at rest-frame 1280 A of AB^q = 19.10 given 
by Fan et al. (2003). The spectral index of the power 
law is a free parameter and the fitting produces a value of 
a = 0.2 (defining /„ oc v~ a , where f v is the flux-density at 
frequency v). This continuum slope is typical of quasars 
(e.g. Vanden Berk et al. 2001). The goodness of fit for the 
Fe II -fcontinuum fit is determined from a chi-square test 
using the noise array from the observation. 

Once the minimum chi-square fit has been determined 
the best-fitting combination of continuum and Fe n emis- 
sion is then subtracted from the spectrum, leaving the iso- 
lated Mgll emission line in the region 2700 A - 2900 A. 
The Mgll emission line itself is then modeled using 
two identical gaussians which represent the line doublet 
(AA2796,2802). During the fitting of the line profile, the 
amplitude, FWHM and central wavelength of the gaus- 
sians are treated as free parameters, while the doublet sep- 
aration is held fixed at its laboratory value. The best fit to 
the Mgll emission was determined with a chi-square test. 
Although more complex dcconvolutions of the line pro- 
file are obviously possible, due to the signal-to-noise con- 
straints of the spectrum the above procedure was adopted 
because it was able to effectively describe the Mgn line 
profile with the minimum level of complexity. 

We also performed an alternative fitting procedure in 
which the continuum was fixed at two points with low 
expected Fe n emission and then the Mg II and Fe II fea- 
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tures fit simultaneously. The results from this indepen- 
dent fitting method are consistent with those using the 
former method. The ratio of Fen/Mgn emission lines 
in high redshift quasars is a powerful diagnostic of early 
element production inside supernovae (Yoshii, Tsujimoto 
& Kawara 1998). Unfortunately we do not have sufficient 
signal to noise to place meaningful constraints on this ratio 
using either of our two fitting methods. 

The best fit Mgn line doublet in SDSS Jl 148+5251 is at 
a redshift of z = 6.41 ± 0.01. The deconvolved FWHM of 
each component of the doublet is 6000 igoo° kms" 1 . The 
rest-frame equivalent width of the Mgn line is 35 A. The 
typical equivalent width of Mgn lines in lower redshift 
SDSS quasars is 32 A (Vanden Berk et al. 2001). The fact 
that the equivalent width of Mgn in SDSS J1148+5251 is 
similar to that of other quasars indicates that the UV con- 
tinuum emission from this quasar is not strongly beamed. 
Fan et al. (2003) noted that the Lya equivalent width 
of this and other z ~ 6 quasars are smaller than lower 
redshift quasars. Pentericci et al. (2002) observed the 
Civ line in the z = 6.28 quasar SDSS J1030+0524 and 
also found an equivalent width typical of lower redshift 
quasars. These results suggest that the low Lya equiva- 
lent widths are primarily due to Ly a absorption and not 
a fundamental difference between the emission properties 
of these quasars. 

3. AN ACCURATE REDSHIFT DETERMINATION 

The redshift of z = 6.43 ± 0.05 for this quasar given in 
Fan et al. (2003) is based on the position of the heavily 
absorbed Ly a line and the onset of the Lyman break. The 
lack of distinguishable Nv emission or any other lines in 
the optical spectrum mean that a more accurate redshift 
could not be obtained. In our spectrum we are able to 
measure the location of the Mg n line and hence obtain a 
more accurate redshift for this quasar. We find a redshift 
of z = 6.41 ±0.01. The Mgn line is particularly useful for 
redshift estimation compared to, for example, C iv which 
could be observed in the J-band. The low-ionization Mgn 
line is usually found close to the systemic redshift, whereas 
the high-ionization C IV line often shows blueshifts of sev- 
eral thousand kms -1 (Richards et al. 2002). Another 
important fact is that there is no systematic velocity off- 
set for the Lya and Mgn lines in a large sample of SDSS 
quasars (Vanden Berk et al. 2001), so we can predict the 
intrinsic, unabsorbed Lya line peak. 

Accurate redshifts for the most distant quasars are im- 
portant for several reasons. The rest-frame UV spectra 
of z ~ 6 quasars show a deficit of flux shortward of Ly a 
due to absorption by neutral hydrogen in the intergalactic 
medium (Fan et al. 2001; Becker et al. 2001). These ob- 
servations have led to the possibility that we are now wit- 
nessing the epoch of reionization. Whilst accurate quasar 
redshifts are not essential for determining the ionization 
state of the universe as a function of redshift, they are 
very important when dealing with the state of the IGM 
close to the quasar. 

Assuming that the IGM in general is neutral at the red- 
shift of the quasar, the strong ionizing field of the quasar 
will ionize a local region around the quasar allowing Ly a 
photons to be transmitted. This proximity effect has been 
observed as a lack of absorption lines in the spectra of 



quasars close to the peak of Lya (Bajtlik, Duncan & Os- 
triker 1988). Cen & Haiman (2000) showed that the size 
of this ionized region and the transmitted flux blucward 
of the peak of Lya depend upon the number of ioniz- 
ing photons which have been emitted from the quasar, i.e. 
the product of its lifetime and luminosity. Haiman & Cen 
(2002) apply this to the spectrum of the z — 6.28 quasar 
SDSS J1030+0524 to show that it is not strongly gravi- 
tationally lensed and to derive a lower limit on how long 
the quasar has been active. Such an analysis depends crit- 
ically on knowing where the peak of the Ly a line should 
be in the spectrum. 

The spectrum of SDSS J1148+5251 presented by Fan 
et al. (2003) also shows zero flux at Lya redshifts of 5.7 
to 6.33, suggestive of a neutral universe. However, there 
is some flux apparent on the blue wing of Lya which is 
likely due to the proximity effect. The marked location of 
Lya at z = 6.43 given by Fan et al. is at the peak of the 
Lya emission. However, adopting the new redshift for the 
quasar of z = 6.41 shifts the location of Lya to 9010 A, 
which is co-incident with a steep decline in flux to less 
than half of the peak value. Hence the revised redshift 
indicates a substantially more neutral IGM close to the 
quasar. The amount and extent of transmitted flux on the 
blue side of Lya is similar to that in SDSS J1030+0524 
and the conclusions about this quasar made by Haiman 
& Cen (2002) and Pentericci et al. (2002) probably apply 
equally to SDSS J1148+5251. 

The evolutionary state of the host galaxies of high red- 
shift quasars can be probed by (sub)millimeter obser- 
vations of dust and molecular gas. Particularly impor- 
tant tracers are the emission lines from transitions of CO 
molecules and carbon atoms (Blain et al. 2000). However 
current millimeter receivers have a fairly small bandwidth 
corresponding to ~ 2000 kms -1 . Hence systemic redshifts 
must be known to an accuracy of ~ 0.01 to ensure that 
the emission line will fall within the frequency range of 
the observation. Our observations of SDSS J1148+5251 
provide a redshift of sufficient accuracy for millimeter line 
searches. 

4. ESTIMATING THE BLACK-HOLE MASS 

Reverberation mapping studies of low redshift quasars 
have been successfully employed to estimate the masses of 
their black holes (Netzer & Peterson 1997). This method 
determines the radial distance of the broad line emit- 
ting region (BLR) from the nucleus (-Rblr), by observ- 
ing the time lag (~ 1 year) between variations in the 
UV/optical continuum and the line strength. Combin- 
ing this distance with the velocity of the emission line gas 
(Ublr, which is closely related to the observed FWHM 
of the line) gives a virial estimate of the black-hole mass 
(Mbh = G _1 i?BLR^B LR ). This estimate assumes that the 
dynamics of the gas are dominated by gravitational forces 
(see Peterson & Wandel 2000 for a discussion). These re- 
verberation studies have revealed a correlation between the 
continuum luminosity and the radius of the BLR which can 
be used to estimate black-hole masses when only the line 
FWHM and luminosity are available (Kaspi et al. 2000). 

Traditionally, the H f3 emission line has been used for 
this work, since it is readily observable in the optical spec- 
tra of low redshift quasars where the reverberation map- 
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ping calibration of luminosity and BLR has been derived. 
Recently, McLure & Jarvis (2002) have cross-calibrated 
this relationship from the H /3 line to the Mg II line using 
a sample of AGN with reverberation mapping measures of 
-Rblr, so that black- hole mass estimates can be derived to 
higher rcdshifts. They note that the Mgil line (which has 
an ionization potential close to that of H/3) is a much bet- 
ter choice than higher ionization lines like C iv whose dy- 
namics may be strongly influenced by outflows. McLure & 
Jarvis derived a calibration between -Rblr and £3000 (the 
monochromatic luminosity at 3000 A) and performed a fit 
to the H/3 and Mgn FWHMs, finding a linear relation- 
ship with virtually identical normalization. The resulting 
black-hole mass estimator has the form 

Mbh _ _ _ 7 / AL3000 \ 47 f FWHM(Mgll) ^ 2 
M - 6 - 67 {w^Wj { kms- 1 ) ■ [ ) 

Our spectrum covers the rest wavelength region 3000 A 
and we have determined the Mg 11 FWHM. Therefore we 
can apply this formula to derive an estimate of the black- 
hole mass in SDSS J1148+5251. Using XL 3Q0Q = 6.2 x 10 39 
Wand FWHM(Mgn) = 6000 kms^ 1 we get M bh = 3xl0 9 
M©. Assuming that z ~ 6 quasars arc not fundamen- 
tally different to low redshift quasars (evidence support- 
ing this is that their spectra are similar), the accuracy of 
this black-hole mass estimator with respect to those from 
reverberation mapping is a factor of 2.5 (la). 

This is an extremely massive black hole to be in exis- 
tence such a short time (840 Myr) after the beginning of 
the universe. Assuming that this black hole is residing in a 
proportionately massive halo (~ 10 13 M©) then standard 
theories of structure formation predict that this is one of 
the rarest, most massive halos to have collapsed by this 
time (e.g. Haiman & Loeb 2001; Fan ct al. 2001). There 
are two caveats to the derived black-hole mass. Since the 
estimator depends on L3000 , if we have over-estimated the 
intrinsic, isotropic luminosity we will have over-estimated 
the black hole mass. If the luminosity function of z ~ 6 
quasars has a steep slope then there are expectations that 
a significant fraction of the SDSS z ~ 6 quasars will have 
their luminosities enhanced by cither gravitational lcnsing 
or beaming (Turner 1991; Wyithe & Loeb 2002a). Fan 
et al. (2003) discussed -fC-band images of the three new 
z ~ 6 quasars and found that none of them show evi- 
dence for gravitational lensing. Our spectrum found that 
the equivalent width of the Mg II line is typical of quasars 
which shows that the continuum luminosity is not strongly 
enhanced by beaming. Therefore, the available evidence 



suggests that there is minimal amplification of the intrinsic 
quasar luminosity. 

It is interesting to consider the maximum luminosity 
that can be attained by an accreting black hole of this 
mass. This is usually considered to be the Eddington limit 
at which the outward radiation pressure equals the inward 
gravitational attraction. The Eddington limit for a black 
hole of mass 3 x 10 9 M© is 4 x 10 40 W. Using a bolo- 
mctric correction of 7 to get the bolomctric luminosity 
from AL3000 (Wandcl, Peterson & Malkan 1999) we de- 
termine a bolometric luminosity for SDSS J1148+5251 of 
L = 4 x 10 40 W. Therefore we find that the black hole in 
this quasar is radiating the maximal amount possible. This 
supports the often made assumption when relating quasar 
luminosities to black hole and halo masses that the high- 
est redshift quasars are accreting at the Eddington limit 
(e.g. Fan ct al. 2001; Wyithe & Loeb 2002b). It is not 
too surprising to find that SDSS J 1148+5251 has the min- 
imum size black hole given its luminosity. As mentioned 
above, the short cosmic time available for the growth of su- 
permassive black holes by this redshift means black holes 
with such masses will only have arisen within rare col- 
lapsed peaks. In such conditions, there is likely to be a 
plentiful gas supply to fuel the quasar and the black hole 
is expected to be accreting close to its Eddington limit. 

5. CONCLUSIONS 

We have presented H and K spectra of the quasar SDSS 
J 1148+5251 which currently has the highest known red- 
shift. We clearly detect the broad Mgn emission line. Fit- 
ting the location of this emission line we derive a redshift 
of z = 6.41 + 0.01. The equivalent width of this line is sim- 
ilar to that of lower redshift SDSS quasars, suggesting that 
the optical continuum is not strongly beamed. We obtain 
a virial estimate of the black-hole mass of Mbh = 3 x 10 9 
M© . The quasar is radiating at the Eddington luminosity, 
in line with the expectation that it is easier to fuel a super- 
massive black hole than create one at such an early epoch. 
Similar observations of other z > 6 quasars would reveal 
if they are all radiating at the Eddington luminosity. 
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